Two-mode squeezed vacuum state coupled to the common thermal reservoir 
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Entangled states play a crucial role in quantum information protocols, thus the dynamical behavior 
of entanglement is of a great importance. In this paper we consider a two-mode squeezed vacuum 
state coupled to one thermal reservoir as a model of an entangled state embedded in an environment. 
As a criterion for entanglement we use a continuous-variable equivalent of the Peres-Horodecki 
criterion, namely the Simón criterion. To quantify entanglement we use the logarithmic negativity. 
We derive a condition, which assures that the state remains entangled in spite of the interaction 
with the reservoir. Moreover for the case of interaction with vacuum as an environment we show 
that a state of interest after infinitely long interaction is not only entangled, but also pure. For 
comparison we also consider a model in which each of both modes is coupled to its own reservoir. 

PACS numbers: 03.65.Ud, 03.67.-a, 03.67.Mn, 42.50.Dv 



I. INTRODUCTION 

A number of possible practical applications of entan- 
gled states has been proposed, including the quantum 
computation [jj and the quantum teleportation 0, 
Thus it is very important to know, how the interaction 
with an environment influences the dynamical behavior 
of entanglement. In many experiments as a source of en- 
tangled pairs serves the process of the parametric down 
conversion in the non- linear crystals 0, , in which two- 
mode squeezed vacuum states are produced. Hence, here 
we consider the time evolution of the two-mode squeezed 
vacuum, as an example of continuous variable entangled 
state, embedded in a thermal environment. Moreover, 
let us assume that, for some reason, e.g. storage, both 
modes are captured in a resonant cavity, i.e. they are 
coupled to a common reservoir. Additionally we pre- 
sume that they have the same frequency. There have 
been intensive research in this direction, but Jeong et 
al. @, Kim and Lee Q, Scheel et al. Q, Hiroshima Q 
and Paris assumed that the two-mode squeezed vac- 
uum interacts with the environment in a diffcrent way. In 
their model both modes are spatially separated and thus 
each of them is independently coupled to its own reser- 
voir. Both modes are taken to have the same frequency 
and the same coupling constant. This model is justified, 
while the two-mode squeezed vacuum can be used as a 
quantum channel between sender and receiver stations in 
the quantum teleportation Q • There is a possible exper- 
imental investigation of degradation of entanglement duc 
to the interaction with the environment . 



II. MODEL AND DESCRIPTION 

In this paper we consider a bipartite system, namely, 
two modes of the electromagnètic field. Each of those 
modes has the Hilbert space equivalent to the Hilbert 
space of the harmònic oscillator. The Hilbert space of 



the system of the interest is a tensor product of sub- 
system's Hilbert spaces, TL = TL\ <g> TÍ2- Each mode is 
described by dimensionless quadrature operators defined 
as Xj = (àj + àj)/V2 and pj = — àj)/y/2 (j = 1,2). 
Those operators obey a Standard commutation relation 



Pk] = iSjk, which is equivalent to 



Ojk- 



Bipartite system in question is embedded in the ther- 
mal environment, which we imitate by the set of har- 
mònic oscillators. We model the coupling of both modes 
to the common reservoir by the interaction Hamiltonian 
of the form: 



ti. 



inl 



b k ). 



(1) 



k=l 



k=l 



While in the case when each mode interacts with its own 
reservoir, the interaction Hamiltonian is: 
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In both cases r/ is a coupling constant. 

In order to study the evolution of our system we use 
the Wigner distribution function defined as [Ï2| : 



(3) 



W(x,p) = tt~ 2 / dV (x - x' | p | x + x') 

x exp(2ix' • p), 



where p is density operator , x = (xi , X2) and p — (pi , ^2 ) • 
For the sake of convenience we arrange the 
quadrature operators and the phase space variables 
into four dimensional vectors X = (x\,pi, X2,p>2) and 
X = (x-i , pi , 0:2,^2)1 and define a quadrature variance ma- 
trix jl3j: 

V aP = ({AX a , AXfj}) = Tr({AX a ,AX p }p) 

= J d i XAX a AXpW(X), (4) 
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where {AX a ,AX } = (AX a AX + AX f3 AX a )/2. Of 
course AX a — X a — (X a ) and (X a ) = Tr{X a p). Sim- 
ilarly, AX a — X a — (X a ), where (X a ) is calculated as 
an average of X a with respect to the Wigner distribution 
W(X.). Let us stress that {X a ) equals to (X a ). 

A (zero-mean) Gaussian state is described by the 
Wigner distribution in a form: 

W(X) = (2^)- 2 (dctV)- 1 / 2 exp(-ÍxV- 1 X T ), (5) 

where V is the variance matrix defined in eq. 10} . 

To verify whether our system is in an entangled state 
we use the generalization of the Peres-Horodecki cri- 
terio n IÏ4L Il5| for continuous variables, derived by Si- 
món [Ï6j- In terms of the quadrature variance matrix, 
eq. (0J), rewritten as follows: 



A C 
C T B 



it takes a form of an inequality |16| : 



detAdetB + |cfeíC|) 2 

-Tr [AJCJBJC T J] - ^(detA + detB) > 0, 
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where matrix J is: 



1 

-1 
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Inequality is a necessary condition on the variance 
matrix of a separable bipartite state 16]. Moreover, Si- 
món |16J has shown that this is also a sufficient condition 
for separability for all bipartite Gaussian states composed 
of two modes, each mode held by each party ^ï}- Thus 
for Gaussian states considered in this paper, when the in- 
equality Q) is violated the state is entangled, otherwise 
it is separable. 

Let us note, that there is possible an equivalent ap- 
proach to the separability properties of quantum two- 
party Gaussian states in the framework of the operator 
formalism for the density operator, as it has been shown 
by Englert and Wódkiewicz |Ï8| . 

To compute the degree of entanglement we use the log- 
arithmic negativity |19| : 



E(p) = log 2 \\ï 



(9) 



where ||/3 Ta ||i denotes the trace norm of p TA , and p TA is 
a partial transpose of the bipartite state p with respect 
to the subsystem A. This measure can be regarded as a 
quantitative version of the Peres-Horodecki criterion 0, 
Ha ] , because it quantifies the degree to which p TA fails to 
be positive [Ï9| . 

For states described by the Wigner function of the 
form JSJ, computing of the logarithmic negativity @ 
simplifies to process known as sympletic diagonalization 



of the variance matrix (0J, as has been shown in [Ï9| . 
i.e. if (Ai, A2) is the sympletic spectrum of the variance 
matrix then the logarithmic negativity equals: 



E{p) = F{\{] 
where F{\) (i = 1,2) is: 
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Evolution of two modes of the electromagnètic field 
coupled to the common thermal reservoir might be de- 
scribed by the diffusion equation, i.e. the Fokker-Planck 
equation. Namely, for the interaction Hamiltonian of the 
form it reads (in the interaction picture): 

d t W(X,t) - Yl [\(P Xj x u + d PjPk ) 

j,k=l,2 

+ ï^±(dl. Xk +dl. Pk )]w(x,t). 



(12) 



where 7 is a coupling constant and N is a mean ther- 
mal photon number. Here it is worthwhile to notice 
that in the case, when one couples each mode to its 
own reservoir, i.e. the interaction Hamiltonian is of the 
form J21 , the relevant Fokker-Planck equation is obtained 
from eq. i|12[l by substitution j = k. 



III. RESULTS 

The process of the parametric down conversion in the 
non-linear crystals is used as a source of entangled pairs 
in many experiments 0.0]- In such process two-mode 
squeezed vacuum states are produced. Thus as an initial 
state of our system we consider the two-mode squeezed 
vacuum, which is described by the Wigner function of 
the form JSJ with the variance matrix as follows: 



V 
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Cl Tl\ 
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where n, and c, are: 



rii = cosh(2r), 
c 2 = —ci = sinh(2r) 



(14) 



where r is a squeezing parameter (for r =/= this state is 
entangled). Without loss of generality we have assumed 
that the squeezing parameter, r, is real. 

While our system evolves, the general form of the vari- 
ance matrix, eq. I|13|) . remains constant, but elements rii 
and Ci are changing. Thus we may rewrite Simón crite- 
rion, eq. 10, in a form: 



1 

16 



[(ni-|ci|)(n 2 -|c 2 |)-l][(ni + |ci|)(n 2 + |c 2 |)-l] > 0. 

(15) 
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Simón criterion, eq. Q, is a yes-no test for entangle- 
ment and thus only the sign of the left hand side of the 
inequality 115(1 is important. Therefore we omit an ex- 
pression, which does not change the sign of the inequal- 
ity (|15f) . Namely, because |cj| is always positive, unless it 
is equal to zero and n\n2 is always greater than one, what 
is a consequence of the uncertainty principle, it is more 
convenient to check the sign of the following inequality, 
instead of using the inequality (|15|) , Q : 



(n x - \ci\)(n 2 - |c 2 |) - 1 > 



(16) 



When it is satisfied then state is separable, otherwise it 
is entangled. 

Solving eq. I(12() with the initial variance matrix in the 
form of eq. 1(13(1 with rii and Ci given in eq. ((14(1 leads to 
time-dependent form of rii and c,-: 



ni = i[2cosh(2r) + (N- e~ 2r )T], 
n 2 = i[2cosh(2r) + (N — e 2r )r], 
ci = \[{N -e- 2r )r- 2sinh(2r)], 
c 2 = è i( N - e 2r ) T + 2 sinh(2r)] , 



(17) 



where r = 1 — e~ 27 *, iV = 27V+ 1 and 7 denotes coupling. 
Then substituting 1(1711 into the inequality 1(161) we find 
that if the initial state is sufficiently squeezed, i.e.: 



|r| > ~ln(2iV+l), 



(18) 



it will remain entangled forever in spite of interaction. 
Otherwise, the state will disentangle after time: 



1 , Í2N + 1 - e- 2 l r l 
t = — ln 



2 7 V 2iV + 1 - e 2 M 



(19) 



Thus basing on the condition ((18|) it is possible to 
choose such initial state of the two-mode squeezed vac- 
uum, for which the entanglement will survive the inter- 
action with the common reservoir. 

Inequality 1(18(1 seems to have the same form as a condi- 
tion for non-classicality of a squeezed thermal state [2CÍj| . 
In the latter case, however, N is a mean thermal photon 
number and r is a squeezing parameter, both correspond- 
ing to the thermal squeezed state. Here N corresponds 
to a number of thermal photon in the environment and 
r is a squeezing parameter of a state embedded in this 
environment. 

For a comparison, if one takes into consideration the 
coupling to two independent thermal reservoirs , it will 
lead to the same general form of the quadrature variance 
matrix, eq. I(13|) . and the separability condition, eq. Q16II . 
but coefhcients n.$ and Cj would be different, namely kJ, 

"ET 
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where, as before, r is a squeezing parameter, N = 2N+ 1 
and 7 denotes coupling. In such a case, every state will 
disentangle after the time: 



1 

7 



■ln 1 



1 



-2|r| 
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unless it would interact with vacuum, as it has been 
shown by Duan et al. , Scheel et al. y| and Paris • 
Such a situation corresponds to the limit N — > 0, what 
implics t — > 00. 

For the variance matrix of the form 1(13(1 it is easy to 
find the sympletic spectrum and the logarithmic ncgativ- 
ity. The sympletic spectrum takes form: 



1 



Ai = -VuH c i)( n 2 + c 2 ), 



A 2 = gVA^ï + c i)( n 2 - c 2 ) 



(22) 



Using above mentioned sympletic spectrum, eq. 1(22(1 . and 
conditions for function F(Xi), eq. 1(11(1 . together with 
eq. ((17(1 it is possible to derive inequality 1(18(1 and ex- 
pression for the disentanglement time in the one reser- 
voir case, eq. l(TÜ|) . Analogously it is possible to derive 
formula for the disentanglement time in the case of the 
interaction with two reservoirs, eq. (|2T1) . 

Especially, sympletic spectrum, eq. ((22(1 . holds for ini- 
tial two-mode squeezed vacuum state with rii and Cj given 
by eq. ((141) , and leads to conclusion that the logarithmic 
negativity for such state is simply propotional to the ab- 
solute value of the squeezing parameter, r: 



(23) 



Moreover, having time-dependent form of rii and a, 
eq. I(17|) . allows us to write down the analytic expression 
for the asymptotic value of the logarithmic negativity at 
t — * 00: 



E( P ) 



ln.2' 



1 



Iog a (2JV + l), 



(24) 



which is vàlid for initial states fulfilling condition l(18() . In 
case, when \r\ = | ln(2jV+l), lim^oo E(p) = (compare 
with event of the interaction of two-mode squeezed state 
with two separate reservoirs in a vacuum state). 

Figs. H an d HI show changes of the logarithmic neg- 
ativity, E(p), during the interaction with the common 
reservoir for different initial vàlues of the squeezing pa- 
rameter, r, and the mean thermal photon number, TV. 
Logarithmic negativity, E(p), is plotted as a function of 
the dimensionless rescaled time r = 1 — exp(— 27Í) which 
equals to zero for t = and to one for t — 00. Note that 
the coupling constant, 7, changes only the time scale and 
has no effect on the logarithmic negativity itself. From 
Fig. ^ h is eas Y t° notice that states satisfying condi- 
tion l(18|) (i.e. for r = 0.5, r = 1 and r = 2) remain 
entangled forever in spite of interaction. 
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FIG. 1: Logarithmic negativity, E(p), for the two-mode 
squeezed state as function of a dimensionless time, r = 
1 — exp(—2 r yt), which is for t = and 1 for í = oo, in the 
common reservoir model, for N = 0.5 and r = (0, 0.1, 0.5, 1, 2) 
(from the bottom to the top). 



FIG. 3: Logarithmic negativity, E(p), for the two-mode 
squeezed state as function of a dimensionless time, r = 
1 — exp(—~/t), which is for t = and 1 for t = oo, in the two 
reservoirs model, for N = 0.5 and r — (0,0.1,0.5, 1,2) (from 
the bottom to the top). 



E(p) 




FIG. 2: Logarithmic negativity, E(p), for the two-mode 
squeezed state as function of a dimensionless time, r = 
1 — exp(—2 r yt), which is for t = and 1 for í = oo, in the 
common reservoir model, for r — 1 and N — (1, 2, 6, exp(2), 9) 
(from the top to the bottom), N = 2N + 1. 



FIG. 4: Logarithmic negativity, E(p), for the two-mode 
squeezed state as function of a dimensionless time, r = 
1 — exp(-'yt), which is for t = and 1 for í = oo, in the two 
reservoirs model, for r = 1 and N = (1, 2, 6, ea;p(2), 9) (from 
the top to the bottom), N = 2N + 1. 



Figs. and 01 show changes of the logarithmic nega- 
tivity, E(p), during the interaction with two separated 
reservoirs for different vàlues of the squeezing parameter, 
r, and the mean thermal photon number, N. Again the 
logarithmic negativity, E(p), is plotted as a function of 
the dimensionless time r = 1 — exp(— jt) which equals to 
zero for t — and to one for t = oo. From Fig. [J|it is 
easy to notice that only in the case of interaction with 
vacuum, namely N = 1 states remain entangled. 

Let us stress that Figs.n a nd|31(similariy Figs.|21and0} 
are plotted for the same initial vàlues of the squeezing 
parameter, r and the mean thermal photon number, N, 
but for different models of interaction. 

It is worth noticing that in a case, when a pure two- 
mode squeezed vacuum state interacts with a common 
reservoir initially in vacuum state (N — 0) the result- 
ing state, after infinitely long interaction, is a pure, but 



it is not a vacuum. Final state is squeezed, and thus 
entangled. 

That is the case because of properties of the interaction 
Hamiltonian Q , which show up after non- local transfor- 
mation of variables, such as: 



(25) 



Were new variables are just sums and differences of the 
old ones. 

Our initial state is pure and rewriting it in new basis 
does not change its purity, while the interaction Hamilto- 
nian Q appears to couple to the reservoir only mode S 
described by sums, {xs,ps)> leaving mode D described by 
differences, (xd,Pd), unaffected [22J. Thus, only mode S 
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FIG. 5: Purity, Tr[p 2 ], for the two-mode squeezed state as 
function of jt, in the common reservoir model, for N = 0.5 
and different initial squeezing. 




FIG. 6: Purity, Tr\p 2 \, for the two-mode squeezed state as 
function of yt, in the common reservoir model, for r = 1 and 
differen mean thermal photon number, N = 2N + 1. 



is not in vacuum, the final state will be mixed. Here 
mode D enables to span some kind of a decoherence-free 
subspace p3|. 

Of course in the case of interaction with two indepen- 
dent reservoirs, the interaction Hamiltonian @ do not 
exhibit such property under change of variables (|25|) . 

For the sake of clarity we plot the purity, Tr[p 2 ], of 
the state of interest as a function of time for different 
initial vàlues of the squeezing parameter, r, and the mean 
thermal photon number, N (see Figs. |S] and [BJ. For a 
Gaussian state described by the Wigner function with 
the variance matrix of the form (|13[) purity equals: 



Tr[p 2 } = 



1 



v / (n 2_ c 2 )(n 2_ c 2 ) 



(26) 



IV. SUMMARY 

Summarizing, for the coupling to the common reser- 
voir, we noticed that whether the state is entangled or 
not depends mostly on the initial degree of squeezing. 
If the state is initially sufficicntly squeezed, i.e. accord- 
ing to condition i|18fl . it will always remain entangled, 
independently of the strength of the coupling to the en- 
vironmcnt. Otherwise, the state initially entangled will 
become separable after the time given in the eq. I|19l) . It 
is clearly different result from that obtained for the model 
in which each of two entangled modes interacts with its 
own reservoir. In such a case only the interaction with 
vacuum fluctuations does not lead to the disentanglement 
|2Ïj . Moreover, we have shown that in a case of interac- 
tion with a common reservoir being in a vacuum state, 
the state of interest, after infinitcly long interaction is 
not only still entangled, but also it is pure. 



interacts with thermal reservoir and mode D undergoes 
free evolution. As a consequence, mode S first experi- 
ences decoherence and then dissipation to the state of the 
reservoir (see Figs. [5] and [ü]- decoherence for jt = — 0.5 
and dissipation for yt > 0.5). Therefore, if the vacuum 
is the environment of interest, mode S becomes vacuum 
after some time. Clearly both modes are pure, mode S 
because it is vacuum and mode D because it had evolved 
freely. Hence the final state is also pure, but it is still 
squeezed, because some squeezing "survived" in mode 
D, and thus it is still entangled. If the reservoir initially 
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